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Abstract 

The Gross-Pitaevskii (GP) equation admits a two-dimensional solitary wave 
solution representing two mutually self-propelled, anti-parallel straight line 
vortices. The complete sequence of such solitary wave solutions has been com- 
puted by Jones and Roberts (J. Phys. A, 15, 2599, 1982). These solutions 
are unstable with respect to three-dimensional perturbations (the Crow in- 
stability). The most unstable mode has a wavelength along the direction of 
the vortices of the same order as their separation. The growth rate associ- 
ated with this mode is evaluated here, and it is found to increase very rapidly 
with decreasing separation. It is shown, through numerical integrations of the 
GP equation that, as the perturbations grow to finite amplitude, the lines 
reconnect to produce a sequence of almost circular vortex rings. 

1 Introduction 

The experimental realization of Bose condensates in weakly interacting gases moti- 
vated the recent explosion in theoretical studies of its properties using the the Gross- 
Pitaevskii (GP) model which represents the so-called mean-field limit of quantum 
field theories. The same equation has been the subject of extensive studies also in 
the framework of superfluid helium at very low temperature. In this case the GP 
model is assumed to be linked to the condensate fraction of the superfluid. This 
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is the ninth in a series of papers devoted to modehng flows in a Bose Condensate. 
Reference wifl be made to the fourth and fifth papers in the sequence: Jones and 
Roberts (1982) and Jones, Putterman, and Roberts (1986). 

Superfluid turbulence has been the focus of many experimental and theoretical 
studies (see Donnelly, 1991). Superfluid turbulence manifests itself as a tangle of 
quantized vortex lines. The dynamics of the tangle depend crucially on the interac- 
tions of the vortex filaments. These have been studied by using the classical model 
of vortices in an incompressible Euler fluid. This omits, however, two mechanisms 
that are very relevant to the superfluid tangle. 

First, as Vinen (2000) argues, emission of sound by a vortex tangle is very sig- 
niflcant in superfluid turbulence. This process is completely removed by the main 
assumption of classical vortex theory: V-v = 0, where v is the superfluid velocity. 
The dynamics of vortex fllaments in a compressible fluid is not as well understood 
as that for the incompressible case. The scattering of sound by compressible Euler 
fluids has however been the subject of several recent investigations; see, for example. 
Ford and Smith (1999). 

Second, the processes of severance and coalescence of vortex lines are centrally 
important for the study of superfluid turbulence, but these are expressly forbidden 
by the Kelvin-Helmholtz theorem, according to which vortex lines are frozen to an 
Euler fluid and cannot change their topology. In an Euler fluid, the processes have 
been successfully simulated numerically by restoring viscosity. This step is disallowed 
in a superfluid, and the only way to defeat the theorem is through ad hoc procedures. 
For example, it was supposed by Schwarz (1988) that, whenever one vortex fllament 
comes within a distance A of another filament, rcconnection will always occur, and 
that otherwise reconnection will not happen. A precise way of determining A is 
not known, but its value can clearly greatly affect the reconnection rate in a vortex 
tangle. Moreover, the angle at which the vortex fllaments approach one another is 
undoubtedly an important factor in determining whether they reconnect or not; a 
clear set of reconnection rules is lacking. 

The advantage of GP theory in comparison with the classical approach is that 
it gives superfluid vortex lines their own unique core structure. At the same time, 
it provides a mechanism for the severance and coalescence of vortex lines, and in- 
cludes sound propagation, so that the acoustic emission from a vortex tangle can be 
evaluated. Koplik and Levine (1993, 1996) used numerical simulations of the GP 
model to study the reconnection of, and the interaction between, straight-line vor- 
tices and vortex rings. In particular, they witnessed the annihilation of vortex rings 
of similar radii. Recently Leadbeater et al. (2001) elucidated the loss of energy to 
sound emission during vortex ring collisions. Their calculations suggested that the 
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sound emitted during reconnections is a significant decay mechanism for superfluid 
turbulence. 

In tliis paper, we first study the hnear stabihty of a vortex pair. This is a two- 
dimensional (2D) structure that can, in GP theory, be represented by a wavefunction 
ipo{x,y,t) that is independent of the coordinate z and which has two zeros at y = 
±|/t, representing vortices separated by a distance h. The phase of ipo increases 
by 271 round one zero and decreases by 27r round the other corresponding, in the 
hydrodynamic interpretation of ipoy^o a pair of antiparallel vortices (sometimes called 
"point vortices") that move uniformly with speed U in the x— direction as a solitary 
wave, i.e., ipo = ^o(a; — Ut,y), where U is obtained by solving the GP equation in 
2D: 

(Here, and frequently in what follows, we use dimensionless variables such that the 
unit length corresponds to the healing length a, the speed of sound is c = 1/ v^, and 
the density at infinity is poo = 1- Later we shall write ipo = uq + ivo, where uq and 
fo are real.) Solutions of this form were first reported by Jones and Roberts (1982) 
who determined the entire of sequence of such solutions and their associated energy 
per unit length S and momentum per unit length V, both of which decrease to zero 
as /i — 0: 

^ = \ j \VM^dV+\j{l-\M^fdV (2) 
^ = [(^0 - l)V^o - (V^o - l)Vro] dV. (3) 

Multiplying (|^) by xdip^/dx and integrating by parts, Jones et al (1986) showed that 









dx 
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dV. (4) 



They located a critical value /i^ ~ 1.7 of /i, at which the sequence lost or gained 
vorticity. For h < he the sequence has no vorticity, although solitary disturbances 
exist as finite amplitude sound waves in which the two minima of {tpol are no longer 
zero. As — > 0, U approaches the speed of sound c and the acoustic solutions merge 
with the phonon branch of the dispersion curve. 

Jones and Roberts did not examine the stability of their 2D solitary waves. It is 
known that the vortex pair in an incompressible Euler fluid is prone to the so-called 
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"Crow instability" (Crow, 1970). Kuznetsov and Rasmussen (1995) proved that in 
the long- wavelength limit, where k is small compared with h, both the vortex pair and 
solitary acoustic solutions are unstable, but they determined neither the boundaries 
of instability nor the wavelength at which the growth rate is maximal. In this paper 
we first solve the linear stability problem for all h with the particular aim of finding 
the growth rate of the Crow instabilities as a function of the separation h. We study 
the subsequent evolution of the instabilities to finite amplitude by integrating the 
CP equations in 3D. This parallels the corresponding analysis by Moore (1972) for a 
classical fluid, but differs in that healing becomes important as the instability brings 
one vortex core close to the other. Unlike the classical case, reconnection can, and 
does, occur so that the final result is a sequence of almost circular vortex rings. 



2 Linear stability of the vortex pair 

We return to the CP equation in the reference frame moving with the vortex pair: 

-2i^ + 2if/g = V^^ + (l-H>. (5) 

We seek solutions of (|^) in the form ip{x, y, z, t) = ipoix, y) + ip{x, y, z, t) where ifj is 
infinitesimal. The resulting linearized CP equation determines the stability of the 
vortex pair. We separate into real and imaginary parts, u and v, and focus on 
separable solutions of the form 

u = u{x, y) exp[crt — ikz] + m*(x, y) exp[crt + ikz]^ (6) 

and similarly for where * stands for complex conjugation; the functions u and v 
are governed by 



VlyU + 2U— + (l - 3mo -vl-k'^ju- 2uqVqv = 2av, (7) 
VlyV - 2U^ +(l-ul- 3v^o -k^)v- 2uoVoU = -2au, (8) 



and, since the perturbation must vanish at great distances from the vortex pair, we 
have 



u 



0, V — > 0, for s = ^/ x"^ + y"^ oo. (9) 



The linear stability problem posed by (0)-(^ has features in common with the 
corresponding classical stability problem analyzed by Crow (1970). In particular. 
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it follows from (0)-(|^) that cr^ is real. This may be demonstrated by introducing 
adjoint variables u and v that obey (j^) and share the same eigenvalue spectrum. We 
multiply (|^ by u, by v, add corresponding sides, integrate over the interior of 
the cylinder s = S, apply the divergence theorem, discarding the resulting surface 
integrals for 5* cxd by an appeal to (|). We then find that u and v must obey 
(0)-(|^), but with 0" replaced by —a. In short, if cr is an eigenvalue of (0)-(|^), so is 
—a. Since all coefficients in (|^) and (||) are real, a and a* are both eigenvalues. Thus 
in all cases cr^ is real. 

The eigenvalues of (0)-® belong to two distinct types of instability, termed the 
symmetric and the antisymmetric modes: 

symmetric: u{—x,y) = —u{x,y), v{—x,y) = v{x,y) 
antisymmetric: u{—x,y) = u{x,y), v{—x,y) = —v{x,y). 

Kuznetsov and Rasmussen (1995) demonstrated that all long wavelength antisym- 
metric modes are stable and all long wavelength symmetric modes are unstable. In 
fact they showed that the dispersion relation for the antisymmetric perturbation is 

= [kUf (l - J-^ < 0, A; 0, (10) 



VU ^ 

and the growth rate of symmetric perturbation is given by 

where £ and V are the energy and momentum per unit length of the vortex pair; see 
(0) and (0). These have been evaluated by Jones and Roberts (1982) and by Jones, 
Putterman and Roberts (1986) for the entire vortex sequence, from the KPl-soliton 
for P — s> to a widely separated pair of vortices for V oo. In the latter case, it 
was found that, in dimensional units, 



.2 



2tx 



m — ha 
a 



(12) 



Vr-^P^h, Ur^^, (13) 

where a is the healing length and a is the vortex core parameter determined numer- 
ically by Pitaevskii (1961) as a ~ 0.38. To compare ([TT| ) with the result obtained by 
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Crow (1970) we rewrite (^) using the cut-off metliod (Saffman, 1992). According to 
tfiis metliod we estimate tlie vortex cut-off parameter S by comparing tlie velocity 
of a ring of radius R given by the cut-off formula 



U 



2n-a6/R 
aS/R 



1 



sin ^9 



WttR 

with the analytical result (Grant and Roberts, 1971) 



K , 4i? 
m 



U 



AnR 



1 8-R 1 I 
In 1 + a]. 

a 



This comparison gives us 



In 25 



a. 



(14) 



(15) 



(16) 



Using (|T6D as a definition of the cut-off parameter for the GP model we can write 

as 



2tt 



In 



h 

2^5 



;i7) 



which together with (|T3]) and (|TT1) implies that, for symmetric modes, 

2 



2 



In 



h 
2^ 



+ 1 



0. 



This establishes they are unstable for all sufficiently large wavelengths. 

It is possible here to compare (|TB|) with the classical theory of Crow (1970), in 
which 6 is the cut-off employed when vorticity is assumed to be confined to filaments. 
Crow assumed the uniform core vortex model, but his derivation is easily adapted 
to a vortex pair with other core structures. He found that, provided ka <^ 1, where 
a is the core radius, 

2 



a 



X 



1 +khKi{kh) - \{khfu;{5) 
1 -khKi{kh) - (khYKoikh) + l{kh)'^uj{6) 



where Kq and Ki are modified Bessel functions, and 

u{6) 



-2 [ (cos u + usinu - 1)^ ^ In(aM) 7 - | + C(5)^ 

JakS 



(19) 



(20) 
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Here 7 ~ 0.577216 • ■ ■ is Euler's constant. When we approximate ( [T^ ) for kh <^ 1, 
we obtain 



1^ 



Now we will address the question of whether the expression ([Tl]) has a more 
general meaning and is valid for the classical core models, so that we can adopt (|18D 
as the general expression for the growth rate of large wavelength perturbations. For 
the uniform core model we relate (0) to the analytical expression for the velocity of 
a vortex ring of radius a: 



U 



K f SR 1 

m 

A-nR V a 4 



This comparison defines the cut-off parameter as 25 
antiparallel uniform core vortices is 



.1/4 



27r 



, h 1 
In - + - 

a 4 



(21) 



The energy of two 
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which, when written using the cut-off parameter, becomes 

.2 
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~27 



, h 1 
In h - 

2a6 2 



(23) 



The momentum and velocity of the vortex pair are given by (|T^). When the right- 
hand side of the expression (|ll|) is evaluated using expressions (|^) and (^) the 
result becomes 



S 



dP/dU 



K 

2Txh 



\ h 1 
In — - + - 

2a5 2 



(24) 



which differs from ([T8|) by |. The nature of this difference together with a brief 
description of the cut-off method are given in the Appendix. 

The linear stability problem (0)-(|^) was solved numerically for various h. The 
region of instability and the maximum growth rate were determined in the fc/i-plane; 
see Figure 1. The stability boundary for the classical hollow core vortices is depicted 
in Figure 1 as well. 
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3 Nonlinear evolution of the instability 



As the unstable perturbation grows in amplitude, it can no longer be described 
by linear equations such as (0)-®- To determine its subsequent evolution it is 
necessary to undertake direct numerical integrations of the GP equation. One can 
then understand the evolution in the following way. As the instability grows, it 
brings some segments of one line into closer proximity with corresponding segments 
of the other line, until the minimum distance between the lines reaches a critical 
value, approximately equal to the critical value [h = he) at which vorticity is lost on 
the solitary wave sequence; see §2. At this moment, the Kelvin-Helmholtz theorem is 
inapplicable and reconnection occurs; curves of zero ip on one line join with the curves 
of zero ip on the other to form closed elongated vortex rings that later relax to become 
approximately circular. Before doing so, each ring oscillates in its fundamental mode, 
being alternatively prolate and oblate; the amplitude of this oscillation diminishes 
as it radiates acoustic waves. 

This scenario is supported by direct numerical simulations, performed with the 
same numerical method as in our previous work (Berloff and Roberts, 2000). In these 
computations we follow the evolution of a vortex pair moving in the x— direction in 
a computational box of dimensions = 60, Dy = 60, Dz = 120. The xy— faces of 
the box are open, to allow sound waves to escape; this is achieved numerically by 
applying the Raymond-Kuo technique (Raymond and Kuo, 1984). The faces z = 
and z = Dz are reflective. To introduce an initial perturbation that does not favor 
any particular wavelength we start with the initial condition 

ipix, y^z,t = 0) = ^poix, y-3)* ipoix, y + 3), (25) 

where 

i)o{x, y) = [l- exp(-0.7r^-^^)] exp(ie) (26) 

is an approximation for the rectilinear vortex and r and 6 are polar coordinates, 
such that X = rcosO and y = rsin^. The wavelength of the instability for h = 6 
was about 30, corresponding to A; ~ 0.2, in good agreement with the result of the 
linear stability analysis, which gave k pa 0.19. This wavelength determines where the 
vortex filaments approach each other and reconnect as vortex rings; see Figure 2. 

The reconnections are accompanied by the emission of sound waves and rarefac- 
tion pulses, resulting in line loss that is approximately one fourth of the total vortex 
line length, so confirming that acoustic losses are significant, and should to be taken 
into account when modeling superfluid turbulence. 
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Figure 2: The isosurface p/ Poo = 0.2 for two anti-parallel vortices initially distance h = 6 
apart that propel each other away from the viewer. An instability develops along the axes 
of these vortex lines, and the lines reconnect to form circular vortex rings. 
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The results of our computations for other h are summarized in Table 1 which 
gives energy, momentum, and wavelength i of the perturbation for the initial field 
and energy, momentum, and radius of the resulting vortex ring and the amount of 
vortex line lost as a percentage of the initial vortex line length. To reduce the time 
taken by an initial perturbation to grow, we took the initial ip to be ipo{x,y+ + 
h/2) * ipo{x,y-h/2) where ipo is given by (|26|) , y± = y ± O.lcosfcz and ipo and k is 
the wavenumber for which the growth rate is a maximum according to the linear 
theory of §2. Notice that if h is small (though larger than the critical distance for 
2D vortex pair) the vortices annihilate each other and the energy is carried away as 
sound waves with the intermediate formation of 3D rarefaction pulses. 

Table 1. 



h 


i 


^init 




'-ring 


V ■ 

' rmg 


Rfing 


% of line lost 


8 


53 


819 


2664 
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14.6 
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201 
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4 Conclusions 

We studied the Crow instability for two antiparallel vortex lines obeying the GP 
equation. Linear stability analysis was used to determine the maximum growth rate 
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of the instability and the region of instabihty. Through numerical simulations of the 
GP equation, we showed that as perturbations grow to finite amplitude the lines 
reconnect to produce a sequence of almost circular vortex rings. We evaluated the 
resulting line loss. 
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Appendix. Classical Crow instability. 

Crow (1970) used the cut-off method to determine the growth rate, cr, of the instabil- 
ity in an incompressible fluid, for all kh and for fca <^ 1, where a is the radius of the 
vortex core. This approximate method is based on the Biot-Savart law determining 
the fluid velocity, v(x), in an incompressible fluid from an assigned vorticity a>(x): 

It is supposed that the vorticity is confined to filaments of infinitesimal cross-section 
so that (p7D reduces to a line integral 

K, f ds' X (x — x(s')) 
^(''' = 4^ / |x-x(.OP ■ (^«) 

where s is arc length on a filament and k is the vorticity contained within it. 

To determine the velocity, U(x), of the filament it is necessary to evaluate (|28| ) for 
each point on the filament x(s), but the resulting integral ( p8|) diverges. In the cut- 
off method the offending segment |s — s'| < a5 of the integral is arbitrarily removed, 
where 6 is the cut-off parameter the value of which depends on the core structure. 
This step is denoted by [6]: 

U(.) = - / ^^;x(>^(^)-^(f)) . (29) 



471- J [5] |x(s) - X( , , 

We give two examples. First for our vortex pair, no divergence arises because ds' 
and x(s) — x(s') are parallel when s' and s are on the same line. Thus 

nh r°° ds K . 
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In the second vortex ring of radius -R(^ a), gives (0). 

An alternative way of defining a cut-off is through the expression for the energy 
of a vortex hne assembly. This is most conveniently expressed as in §153 of Lamb 
(1945) as 

„ off co'(x)'u;(x') , , , , , 

E= / ) ' ^ ' dxrfx, 31 

87ry y |x-x'| ^ ^ 

which, when the vorticity is concentrated into filaments is, 

Stt J[s]Ji5] |x(s)-x(s')| ^ ' 

and [S] signifies that the segment |s' — s| < a6 is removed. Returning to our two 
examples, (|32| ) gives 

Ztt 2ao 



for the vortex pair and 



for the thin-cored ring. 

The cut-offs 6 and 6 must be chosen differently. In order that the Hamiltonian 
relation 

U = dE/dP (35) 

is obeyed by the ring, where P = phiTrP? is its momentum (impulse), it is necessary to 
hold the volume 27r^ i?a^ of the ring constant in the differentiation (p5D; see Roberts 
and Donnelly (1970) and Roberts (1972). This requirement, which implies that 

H5/5) = i, (36) 

is relevant even for the hollow core vortex, since any change in volume would imply 
that pdV work is done on the system at infinity, with a concomitant change in E 



that would cause (|35|) to fail. Relations (|33|) , (|3^) , and (p6D show that 

£ = ^(l„A + n, (37) 



27r V 2a6 2 

a6 



E = -p«:^i?(ln^-- ). (38) 
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Result (^) agrees with known facts for the uniform core (In 25 = i) and hollow core 
(In 25 = |) rings; see Lamb (1945) and Saffman (1992). (The uniform core vortex is 
one in which u/s is constant where ulfj) is the vorticity and 1,^ is the unit vector in 
the direction of increasing 0; see §165 of Lamb, 1945.) Similarly for the vortex pair, 
U = d£/dV, where 

V = pnh (39) 

is the momentum per unit length. 

The expressions of Crow for the growth rates of his antisymmetric and symmetric 
modes of instability reduce, in the limit kh 0, to (0) and (0) above, but not 
with the expression (^) for £. In its place stands 

8=P^LJL + i\. (40) 



27r V 2a5 

It seems to us that this puzzling difference may be connected to the different frame 
of reference used in deriving (^01). Crow used the co-moving frame, in which (^Dp 
translates to the energy £ = £ — VU : 

^=^ln^. (41) 
27r 2a8 ^ ' 



For the reason why this does not coincide with (|37|) we offer the following speculation. 

In the laboratory frame (the frame in which the fluid is at rest at infinity) the 
streamf unction ilj{x,y) for a hollow core vortex is 

V^ = fln^, (42) 

where ri [r2] is a distance from (0, |/i) [from (0, —^h)]; since these distances change 
with time, ip is implicitly a function of t. The streaklines, i.e., the curves parallel 
to the instantaneous direction of the velocity v = V x (iplz) are coaxial circles 
surrounding the cores, the surfaces of which are (for a <^ h) ri = a and r2 = a. The 
energy per unit length is (see §157 of Lamb, 1945) 

£ = hp v'^d^x = lpK{ijjsi - 1ps2) = PKipsi, (43) 



where il^si {'ips2) is the value of on the core surface, ri = a (r2 = a). This correctly 



reduces to (^) for In 25 = | . 
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Consider now the flow as seen in the co-moving frame. This consists of two 
parts: an interior region composed of (non-circular) streamlines surrounding the 
vortices and an exterior region where the steamlines start and finish at oo, where 
V = —Ulx- The two regions are separated by an oblate dividing streamline, x = 
±[tanh(?///i) — y^ — ^h'^Y^'^, on which ip = 0; see the Figure on p. 221 of Lamb (1945). 
Here 

i;=-Uy + ^ = -^(j- + \n'-^^ (44) 

is the streamfunction in the co-moving frame. 

Since the interior fluid is, as seen in the laboratory frame, perpetually carried 
along by the vortex in its motion, it has a special significance. Its energy is [cf. (^3])] 

p«^^, = f!(l„/^-i)=f!|„^=f. (45) 
In \ a 2 / In 2ao 

The same argument applies with minor modifications to vortices with other internal 
structure. 
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